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Abstract. We study the critical energy and magnetization profiles for the Ising 
quantum chain with a marginal extended surface perturbation of the form A/y, 
y being the distance from the surface (Hilhorst- van Leeuwen model). For weak 
local couplings, A < A c , the model displays a continuous surface transition with A- 
dependent exponents, whereas, for A > A c , there is surface order at the bulk critical 
point. If conformal invariance is assumed to hold with such marginal perturbations, 
it predicts conformal profiles with the same scaling form as for the unperturbed 
quantum chain, with marginal surface exponents replacing the unperturbed ones. 
The results of direct analytical and numerical calculations of the profiles confirm the 
validity of the conformal expressions in the regimes of second- and first-order surface 
transitions. 



1. Introduction 

The methods of conformal invariance have been extensively used in two- 
dimensional (2D) critical systems, hrst to study the bulk critical properties [1] but 
also, very soon after, to explore the local critical behaviour at surfaces [2] or near 
line defects [3]. Although translation invariance is broken in the transverse direction 
for such defective systems, conformal invariance is preserved. In the case of the 
line defect, the spectrum-generating algebra was identified as the U(l) Kac-Moody 
algebra [4]. More surprisingly, some aspects of conformal invariance, such as gap- 
exponent relations, seem to hold for marginal extended perturbations, provided the 
shape of the perturbation is properly transformed [5-10]. 

Extended perturbations were hrst considered at the surface of the semi-infinite 
2D Ising model by Hilhorst and van Leeuwen (HvL) [11]. This problem, as well 
as its extensions to linear and radial defects in the bulk, became an active field of 
research during the years that followed [12-20]. Different approaches, either numerical 
or analytical, have been used to elucidate the influence of such perturbations on the 
local critical behaviour (see [21] for a review). 
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Surface extended perturbations correspond to smoothly varying inhomogcneties 
for which the nearest-neighbour coupling, at a distance y from the surface, deviates 
from its bulk value by a quantity Ay~ K . When the decay exponent k is equal to 
the scaling dimension of the couplings, the perturbation is marginal and the local 
critical exponents vary continuously with the perturbation amplitude A [22-24] . For a 
sufficiently strong enhancement of the perturbation amplitude A, the surface remains 
ordered at the bulk critical point. When k is larger than the marginal value, the 
perturbation decreases sufficiently fast and the critical properties remain unchanged. 
On the other hand, for a value of k smaller than the critical one, the system is driven to 
a new surface fixed point, at which power laws are replaced by stretched exponentials 
for A < whereas the surface remains ordered at criticality when A > 0. 

In the present paper, we consider the off-diagonal profiles [25] of the energy and 
magnetization densities for the 2D Ising model with a marginal HvL surface extended 
defect. We work in a strip geometry resulting from a logarithmic conformal mapping 
of the half-space. We treat the problem in the extreme anisotropic limit [26], where the 
2D classical Ising model corresponds to a one-dimensional (ID) Ising quantum chain 
in a transverse field. The energy density profiles are obtained analytically, working in 
the continuum limit, whereas the magnetization profiles are studied numerically. The 
results are compared with the conformal predictions in the different regimes. 

The paper is organised as follows. In section 2, we review briefly the conformal 
results for the surface correlation functions in the presence of a marginal HvL 
perturbation and the resulting off-diagonal profiles in the strip geometry. In section 3, 
we present the HvL model in the extreme anisotropic limit and the diagonalization 
of the associated Hamiltonian in the continuum limit. The off-diagonal profiles are 
calculated in section 4. The results are discussed in section 5. 



2. Marginal extended perturbations and conformal invariance 

Following [10] let us first consider the behaviour, under the infinitesimal special 
conformal transformation z' = z + ez 2 , of a marginal inhomogeneous coupling of 
the form A(p,8) = Ap~ Vx f(9), conjugate to the field x- We use polar coordinates, 
z = pcxp(i9), and y x is the scaling dimension of A. Up to 0(e) we have 

p' = p + ep 2 cos6 6' = 6 + epsin9 (2.1) 

so that the half-plane, with < 8 < ir, is mapped onto itself. The local dilatation 
factor is b(p,9) — |dz'/dz| _1 = 1 — 2epcos8 and leads to the following transformed 
perturbation: 



A(p',6') = (I -2ep cos 6) y xA^ 9 ^ 



pVx 



„ d In f 
1 - ep sin 6 ( y x cot 6 + 



The perturbation is invariant when the coefficient of e vanishes, i.e., when A(p,8) = 
A\ps'm6\~ v * , which is just the form of the marginal HvL perturbation. 

If conformal invariance holds with such a marginal perturbation, the two-point 
correlation functions satisfy the same set of partial differential equations as the 
unperturbed ones [2]. 
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The transformation properties of the two-point functions under the infinitesimal 
special conformal transformation imply the following scaling behaviour at critical- 
ity [2]: 

G^{ Xl -x 2 , yi ,y 2 ) = ( yi y 2 )-**g^) u = ^ ~ ^ + (i/i ~ ^ (2 . 3) 

V1V2 

where x (y) is the direction parallel (perpendicular) to the surface and the bulk 
scaling dimension of the field 4>. Using polar coordinates and taking the limit p 1 ^> p 2 , 
simple scaling considerations lead to 

(\ — x % 
(sin0 lS in0 2 )*J-** (2.4) 

where x s ^ is the surface scaling dimension of <j>. Using the logarithmic conformal 
transformation, w = (L/tt)Id.z, with z = x + iy = pexp(i#) and w = u + iv, 
the half-plane, y > 0, with uniform boundary conditions is mapped onto a strip 
(-co < u < +00, < v < L) with the same boundary conditions on both edges. The 
correlation functions in the two geometries are related under the conformal mapping 
so that 

G^w^) = KzM^G^zJ (2.5) 

where b(z) = |dw/dz| _1 is the local dilatation factor. Thus one obtains the well known 
result for the two-point function in the strip geometry [27] 

G> K -u 2 ,v 1 ,v 2 ) ~(^-J |sin(^— J Sm \~r)\ exp [~L X d u i - u 2)\ ( 2 - 6 ) 
when u 1 3> u 2 . 

This result is expected to apply with a marginal extended perturbation, provided 
the shape of the perturbation is also transformed according to the logarithmic 
mapping. Details will be given in the next section. Then x^ is the continuously 
varying surface scaling dimension of 4> associated with the marginal defect. 

In the strip geometry, the row-to-row transfer operator T of the isotropic classical 
system can be used to expand the correlation function over its eigenstates \n) 
such that T\n) = exp(— E n )\n). This leads to 

G^K -u 2 , Vl , v 2 ) = ^(O|0( Ul )|n}(n|0( U2 )|O}cxp[-( J B„ - £ )K ~ « 2 )1 ( 2 -7) 

n>0 

where |0) is the eigenstate corresponding to the largest eigenvalue of the transfer 
operator and <j> the operator associated with the classical field (j). In the limit u l 3> u 2 , 
the sum is dominated by the eigenstate \<j)) with the largest eigenvalue exp(— E^) and 

a non- vanishing matrix element (0|</>|</>) so that finally one obtains: 

G^C**! - «2, «i, « 2 ) ^ <O|0(wi)|0><^|0(« 2 )|O> exp[-(^ - £? )(«i - « 2 )] ■ ( 2 -8) 

A comparison to equation (2.6) allows us to deduce the off-diagonal profile for <p m 
the strip geometry [25]: 



<t> Q M = mm ~ ( z ) [ sin ( T j] • ( 2 -9) 
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One also recovers the gap-exponent relation [27] 



E = 



(2.10) 



which is valid under this form for the isotropic strip. 

Off-diagonal profiles contain singular contributions only, which facilitates the 
scaling analysis in numerical calculations. Furthermore, an external symmetry- 
breaking field is not needed to obtain a non-vanishing order parameter density. 
These profiles also give an information about the surface critical behaviour since 
equation (2.9) with d<L leads to (j>od( v ) ~ L~ x *. 

One may note that the diagonal profile 



sin ^ 



irv\ 

TJ 



(2.11) 



obtained in [28] for the order parameter, with symmetric fixed boundary conditions, 
is formally recovered by setting = in equation (2.9), i.e., the scaling dimension 
of the order parameter at the extraordinary surface transition. 



3. Diagonalization of the HvL Ising quantum chain in the continuum limit 

Let us start with the classical HvL Ising model on a semi-infinite square lattice 
with nearest-neighbour couplings 

AK 2 (y) = K 2 (y) - K 2 (oo) = Ay~ K K 2 (oo) = K 2 (3.1) 

in the direction perpendicular to the surface and constant couplings, K 1 , in the 
direction parallel to the surface. Hereafter we consider the marginal situation where 
k = 1/v = 1 for the 2D Ising model. Under the logarithmic conformal mapping, the 
marginal perturbation in the half-plane transforms into 



AK 2 (v)=A 




in the strip geometry where < v < L. 

The spectrum of the row-to-row transfer matrix has been obtained analytically [5] 
in the extreme anisotropic limit where K 2 — ► 0, K 1 — > oo while K 2 jK\ = A remains 
constant [26]. Here K\ is the dual coupling such that tanhifj" = exp(— 2K 1 ) and 
self-duality allows to locate the critical point at A c — 1. The amplitude A — > 0, too, 
and thus can be parametrized as A = aK*. 

In the extreme anisotropic limit, the fluctuations are anisotropic with a correlation 
length ratio given by [29] 



C 2 _ cosh(2X 2 ) 



^ cosh(2 J ft' 1 ) 



2K* (3.3) 



Isotropy is restored by rescaling the lattice parameter a 1 in the time direction to 
a 1 = 2K\ measured in units of a 2 in the transverse direction. The row-to-row 
transfer operator, which can be written as T = cxp(— 2K^H), has to be applied 
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l/a 1 times in order to ensure a transfer by one unit of a 2 on the rescaled isotropic 
system. Thus one may identitfy the transfer operator T of the isotropic system to 
exp(— TL) where H is the Hamiltonian of the Ising quantum chain in a transverse field 
such that H\n) = Ejn). 

The critical Hamiltonian reads 



H 



L-l 



.1=1 1=1 

where the cts are the Pauli spin operators, and the couplings A; are given by 
A, = 1 + a 



(3.4) 



' L ■ 1 




— sin 


:f)i 


7T < 





(3.5) 



After a Jordan- Wigner transformation [30], the Hamiltonian (3.4) becomes a 
quadratic form in fermion creation and annihilation operators. It is diagonalized 
through a canonical transformation [31,32] leading to 



^ = E A * (4% -I) • 



(3.6) 



The r/ k s (rj k s) are diagonal fermion creation (anihilation) operators and the A fe s arc 
the energies of the fermionic excitations. They satisfy the following set of equations: 



(A + B)0 fc = A k tp k (A - B)V> fe = A fe fe 



(3.7) 



where <f>k an d ipk are normalized eigenvectors. Working with the xp k s such that 
(A + B)(A - B)tp k = A^fc and changing tp k (l) into (-l) l tp k (l), one is led to the 
eigenvalue problem, 



- \_Mi i) + [\f + i]^ fc (i) - x^ k (i + 1) = a 2 mi) 

with the boundary conditions 

V fe (0) = X^ k (L)-^ k (L + l)=0. 



I = l,L 



(3.8) 



(3.9) 



To study the low-energy spectrum, one may solve the finite-difference equation (3.8), 
with A ; given by (3.5), by working in the continuum limit [5]. To second order in 
1/i, using the continuum variable z = ttI/L, the difference equation is replaced by 
the differential equation 



LA, 



TV 



a cos z 



sin 2 z 



il> k {z) = 



with the boundary conditions 



smz 



o. 



(3.10) 



(3.11) 
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Using the change of function ip k (z) = cos~"(z/2) sin ~ a (z / 2) f k (t) where t = 
sin 2 (z/2), (3.11) becomes a hypergeometric differential equation. Its general solution 
is given by [5] 



f k (t) = C 1 F(a, 6; c; t) + C 2 t L - c F(a - c + 1, b - c + 1; 2 - c; t) (3.12) 



where 



a = -a+- + — — b = -a + - — c=%-a. (3.13) 

2 7T 2 7T 

Due to the boundary conditions in equation (3.11), one has to separately consider 
two different regimes: 

(i) When a < j , the condition at z = imposes C 2 = while the quantization of the 
spectrum, 

A k = j{l- a + k ) k = 0,1,2,... (3.14) 

follows from the condition at z = n. The spectrum exhibits a conformal tower- 
like structure but, to our knowledge, until now no underlying algebra has been 
identified. Using the gap-exponent relation in equation (2.10), one deduces the 
critical dimension of the magnetization from E a — E = A and the critical 
dimension of the energy density from E € — E = A + A 1 . Thus we have 

x s m = \-a x% = 2-2a. (3.15) 

These values are in agreement with the result obtained on a semi-infinite system 
for the magnetic exponent [11] and through finite-size scaling on a strip for both 
exponents [19], in the extreme anisotropic limit. 
Using equations (3.12)-(3.14), one obtains the eigenfunctions 

* (z) ~ 1 sinl '"(^/ 2 ) fc! pG- a >-l- a h COBZ ] (3 16) 
W) - L y 2 cos o (z/2) (3/2 _ a) /k i C0 ^J . (3.16) 

where P^' n \z] is the Jacobi polynomial of degree k [33]. The prefactor, 1/L 1 / 2 , 
follows from the normalization. 

(ii) When a > \, the boundary conditions lead to C 1 = and 

A = O(L- 2a ) A fc = j(i+a + *-l) fe=l,2,3,.... (3.17) 

Again a tower-like structure appears. One has to notice that the differential equa- 
tion (3.10) together with (3.11) gives A = since, in the continuum limit, the 
excitation spectrum is calculated to the first order in 1/L. But, as shown in [5], 
the correct 0{L~ 2a ) behaviour can be extracted from the difference equation (3.8) 
using the trial eigenvector if> (l) ~ tan a (7r//2i) given by (3.10) with A = 0. 
The anomalous behaviour of the first excitation is linked to the persistence of a 
spontaneous surface magnetization at the critical point for a > \ [12,15]. The 
surface order is due to the enhancement of the couplings near the surface which 
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overcompensates the missing bond. The surface displays a first-order transition 
at the bulk critical point with a jump of the surface magnetization to zero when 
A < 1 since, considering the classical picture, the ID surface cannot remain ordered 
when the bulk is disordered. 

From the gap-exponent relation (2.10), the following critical indices are deduced: 



< = o 

The eigenfunctions read 
lf. (*)~-^tan Q (f) 



1 

2 +a - 



(3.18) 



feiW ~ JT72 tan " (I) cosl+2 " Q F \- k + 1 ' k + 2a > \ + a > sin2 (I) 



, (3.19) 



where the hypergeometric function, F, reduces to a polynomial of order k — 1. 
The anomalous normalization of ip (z) is a consequence of the spontaneous surface 
order. 

In analogy with surface transitions in homogeneous systems, (i) is a regime of 
ordinary surface transition. For (ii) we use the name of first-order surface transition. 



4. Off-diagonal conformal profiles 



4-1. Energy- density profiles 

First we consider the energy-density profile for which we derive exact expressions. The 
off-diagonal profile is given by the matrix element e(l) — (0\o~f |e) where |e) = fyj^olO) 
is the lowest excited state with an even number of excitations. Writing a z in terms of 
diagonal fermions, in the continuum limit, one obtains [19] 

e(z) = V' (^)^i(tt -z) + ip (ir - z)^(z) . (4.1) 



(i) At the ordinary surface transition, a < i, using equation (3.16) with k = and 
k = 1, a straightforward calculation leads to the following expression: 

e(z)^y(smz) 1 - 2a (4.2) 

where A is a constant amplitude. This expression agrees with the conformal result 
in equation (2.9) since z = irl/L, the bulk scaling dimension of the energy is x c = 1 
for the Ising model and the surface one is xl = 2 — 2a according to (3.15). 

(ii) At the border line, a. = \, the surface magnetization vanishes logarithmically [12] 
and the first excitation scales normally as 1/L. One may use equation (4.1) 
to calculate the energy profile. In this case both (3.16) and (3.19) give 
the same solution. In particular ip (z) ~ L -1 / 2 tan 1 / 2 (z/2) and ip 1 (z) ~ 
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Figure 1. Rescaled energy-density profile e(z) in the regime of ordinary surface 
transition with a = — ^ for chain sizes L = 2 6 (diamond), 2 7 (triangle), 2 8 (square) 
and 2 9 (circle). The conformal profile (full curve) is shown for comparison. 
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Figure 2. Rescaled energy-density profile e(z) in the regime of first-order surface 
transition with a = 1 for chain sizes L = 2 6 (diamond), L = 2 7 (triangle), 2 8 (square) 
and 2 9 (circle). The result obtained in the continuum approximation (full curve) is 
shown for comparison. 



L 1 / 2 tan 1//2 (2;/2) cos 2 (2;/2). Inserting these eigenvectors into (4.1) one obtains 
a constant profile, 



e(z) = 



D 



(4.3) 
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Figure 3. Rcscalcd energy-density profile ei2(z) in the regime of first-order surface 
transition with a = 1 for chain sizes L = 2 6 (diamond), L = 2 7 (triangle), 2 8 (square) 
and 2 9 (circle). The conformal profile (full curve) is shown for comparison. 



in the continuum limit. This may be traced to the fact that, at a = |, the energy- 
density has the same scaling dimension at the surface as in the bulk. In the energy 
sector the system does not feel anymore the surface and behaves asymptotically 
as a periodic one. The profile (4.3) is recovered using (2.9) with = x c = 1. 
Actually, the numerical data display logarithmic corrections. Good fits of the 
finite-size results were obtained with the functional form e(z) = A Q + A l ln(sin2;) 
where A Q 3> A 1 . 

For a > |, i.e., in the first-order regime, equation (4.1) leads to the following 
expression for the energy density: 



e(z) 



D 



£1/2+0 



cos 



2a+l 



(!H° 2 " +I © 



(4.4) 



Near the surface e(z) scales as L _ ( 1 / 2+a \ i.e., according to (3.18), the gap- 
exponent relation is satisfied but the profile has a functional form which clearly 
disagrees with the conformal expression in (2.9). This is due to the localization near 
z = 7T of the lowest eigenstate ipo( z )i with vanishing eigenvalue in the continuum 
limit. This anomalous behaviour is related to the persistence of surface order 
at criticality. The next two eigenstates are extended and have a normal scaling 
behaviour. Let us consider the corresponding off-diagonal energy profile, 



e i2<» = ^i(2)^ 2 (tt -z)+ip 1 (ir- z)ip 2 {z) . 
Inserting (3.19) with fc — 1,2 yields: 

e 12 (z) = j(sinz) 2a + 1 . 



(4.5) 



(4.6) 
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This is to be compared with (2.9) with x e = 1 and the approriate surface scaling 
dimension, (x*) 12 , which is not given by (3.18) but can be deduced from the gap- 
exponent relation (2.10) and reads: 



«)i 2 = -(A 1 +A 2 ) = 2 a + 2. 



(4.7) 



Thus we find an agreement with the behaviour expected from conformal invariance, 
if we leave out the anomalous first excitation. This point is discussed further in 
the final section where a parallel is made with the analogous situation for the 
homogeneous system with fixed boundary conditions. 

The energy density profiles obtained numerically in the two regimes using 



equation (4.1) for a < | and equations (4.1) and (4.5) for a > 
figures 1-3, respectively. 



i are shown in 



4-2. Magnetization profiles 

The off-diagonal magnetization profile is given by the matrix element m{l) = (0|<t^|o"), 
where \a) = ?7q|0) is the lowest excited state in the odd sector of H with energy 
E a . Rewriting af in terms of diagonal fermions and using Wick's theorem, the local 
magnetization can be expressed as a determinant [34] 



m(l) 



H 



G 



21 



G 1 
G 



Hi G, 



G 



12 



G 



u-i 



with 



(4i 



(4.9) 



where the vectors <p k and ip k are the eigenvectors defined in (3.7). Since erf is non- 
local in terms of the diagonal fermion operators, rjl and n k , the expression of m(l) is 
much more complicated than the corresponding one (4.1) for the energy density and 
we had to study the behaviour of the magnetization profile numerically 

The rescaled profiles shown in figure 4, obtained for different values of a, converge 
with increasing sizes towards the conformal expression which, according to (2.9) with 
i m = | for the bulk magnetic exponent, takes the form 

A 



m(z) 



L 1 / 8 



(sinz) 3 



-1/8 



where z — nl/L. 

For a < i, the critical exponent x s n 



(4.10) 



a has to be used. In particular, the 



value a — | corresponds to equal surface and bulk scaling dimensions, = x m = ^, 
and then (4.10) leads to a constant magnetization profile. The deviations near the 
surfaces in figure 4, which are slowly decreasing with increasing sizes, are due to 
finite-size corrections. 

In the first-order regime, a > i, the critical profile has the same form as for the 
homogeneous system with fixed boundary conditions since x^ n = in both cases. The 
rescaled profile is divergent at the surface in the continuum limit, the numerical data 
obtained on the discrete system display this tendency with increasing sizes. 
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Figure 4. Rescaled magnetization profiles with a = |, a = | and a = — i from top 
to bottom for chain sizes 2 6 (triangle), 2 7 (square) and 2 8 (circle). The full curves 
give the corresponding conformal profiles. 



5. Discussion 

The anomalous behaviour of the energy-density profile in the first-order regime, 
a > |, is linked to the localization, near z = 0, of the state ?p {z). The corresponding 
excitation, A , vanishes as L~ 2a , i.e., faster than higher excitations which display 
the usual L" 1 scaling behaviour. The localization of the lowest eigenstate is itself 
related to the persistence of surface order at criticality. This suggests an analogy 
with the homogeneous system with symmetric fixed boundary conditions for which, 
similarly, the lowest excitation vanishes and corresponds to a mode which is completely 
localized at the surface. The energy spectrum is doubly degenerate and the first 
mode does not contribute to the conformal profile which is given by the matrix 
element (0|cr z |e ++ ), where |e ++ ) = ??2?7i|0) 1S the nrs t non-trivial even state [25]. The 
spectrum and eigenvectors are those of the dual chain with free boundary conditions. 
In the continuum limit, the energy-density profile is given by (2.9) with x e = 1 and 
(x s e ) ++ =L/w(A 1 +A 2 )=2, i.e., 

e ++ {z) = — sinz . (5.1) 
±j 

This conformal profile corresponds to e 12 (z) in the HvL model. The two systems differ 
only by the values of the surface scaling dimensions for the energy: to (2^)++ = 2 
in the model with fixed boundary conditions corresponds x\ = 2 + 2a in the HvL 
model. The difference comes from the shift between the two excitation spectra, 
A fe = (A fe ) ++ + a. Since in both cases the surface perturbation docs not change 
the bulk critical behaviour, the bulk exponents appearing in the profiles are those of 
the homogeneous Ising quantum chain. 

The surface magnetic exponent, x^, vanishes for the two problems due to surface 
order. Thus, if conformal invariance holds as indicated by the numerical data, 
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one expects the two magnetization profiles to differ only through the value of their 
amplitude, which is a-dependent in the HvL model and, asymptotically, we have 

10 ... 

(^2l)++ (6*22)++ ••• (^2Z-l)++ 

... . ( 5 - 2 ) 

(^a)++ (G l2 ) ++ ... (G u _ 1 ) ++ 

where (G i j) ++ is defined as in (4.9) and the index H — h refers to the homogeneous 
chain with fixed boundary conditions. This expression has to be compared to (4.8) 
where 4> (l) = 0(1) and 4> (n) = 0(L~ a ) when z — nn/L = 0(1). All the information 
about the inhomogeneity is then contained in the non-universal amplitude A{a) which 
is increasing with a, like the surface critical magnetization. 

To summarize, we have studied the critical energy and magnetization profiles 
of the HvL model, using the methods of conformal invariance and through direct 
analytical and numerical calculations. We have shown that the scaling form of 
the profiles, resulting from conformal invariance, remains valid for the marginally 
perturbed system. In the regimes of second- and first-order surface transitions, the 
profiles keep the same form as for the unperturbed system with free and fixed boundary 
conditions, respectively, with the HvL values for the surface exponents. 



m(l,a) ~ m ++ (l) = 
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